Size, Shape and Low Energy Electronic Structure of Carbon Nanotubes 
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A theory of the long wavelength low energy electronic structure of graphite-derived nanotubules 
is presented. The propagating 7r electrons are described by wrapping a massless two dimensional 
Dirac Hamiltonian onto a curved surface. The effects of the tubule size, shape and symmetry are 
included through an effective vector potential which we derive for this model. The rich gap structure 
for all straight single wall cylindrical tubes is obtained analytically in this theory, and the effects of 
inhomogeneous shape deformations on nominally metallic armchair tubes are analyzed. 
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Since the discovery of a new family of carbon based 
structures formed by folding graphite sheets into compact 
tube-shaped objects, there has been interest in the elec- 
tronic properties which can be realized with these struc- 
tures Q. It is now understood that these tubes exhibit 
insulating, semimetallic or metallic behavior depending 
on the helicity of the mapping of the graphite sheet onto 
the surface of the tube Discrete microscopic de- 

fects, in the form of disclination pairs, provide an inter- 
face between neighboring straight tubule segments of dif- 
ferent helicities with different electronic gaps, providing 
a novel class of elemental heterojunctions ||. 

In this Letter we investigate the effects of shape fluc- 
tuations on the electronic properties of the carbon nan- 
otubes. We present a new formulation of this problem 
which allows us to study the effects of geometry on the 
quantum dynamics for a tt electron propagating within 
the surface of the wrapped graphite sheet. We show that 
the very rich phenomenology already well established for 
straight single wall cylindrical tubules can be derived di- 
rectly from this geometrical theory We then ex- 
tend the model to consider the effects of inhomogeneous 
deformations in the form of local twists and bends of the 
tubule on the low energy electronic structure. These are 
important low energy structural degrees of freedom of the 
tubules, and indeed one finds that these deformations are 
easily quenched into any three dimensional network com- 
posed of tubules. We show that these shape fluctuations 
also have a very strong effect on the low energy electronic 
and transport properties. 

An isolated two dimensional sheet of graphite is a 
semimetal, with the Fermi energy residing at a critical 
point in the two dimensional tt electron spectrum. The 
Fermi surface is collapsed to a point for this system; there 
are two distinct Fermi points at K (K') points of the zone 
(±47r/3a, 0), where a is the length of the primitive trans- 
lation vector (a = \/?>d where d is the nearest neighbor- 
bond length of the graphite lattice). Expanding the 7r 
electron Hamiltonian around either of these points, and 
linearizing in spatial derivatives, one finds that the low 
energy electronic states are described by a massless two 
dimensional Dirac Hamiltonian, H c g = va ■ p, where p 



denotes a two dimensional momentum in the graphite 
plane, and the ct's are the 2x2 Pauli matrices Q. Here 
the two spin polarizations of the particle refer to the two 
independent basis states (labeling the a and b sublat- 
tices) in the graphite primitive cell. Thus, in addition 
to its physical spin and momentum, the tt electron car- 
ries an internal pseudospin index, labeling the sublattice 
state, and an isospin index, labeling the two independent 
Dirac spectra derived from the K and K' points of the 
zone. The Fermi energy for this system is at E=0. 

To study the electronic behavior of the tubule one 
maps H c g onto a curved surface. The mapping of the 
graphite sheet onto the cylindrical surface can be speci- 
fied by a single superlattice translation vector Th which 
defines an elementary orbit around the waist of the cylin- 
der. In the absence of disclinations the superlattice vec- 
tor Th is an element of the original graphite triangular 
Bravais lattice. This wrapping is conventionally indexed 
by two integers [m, n] such that T% — mT\ + nT 2 with 
Ti = o(l,0) and T 2 = a(l/2, V§/2) |-|. In figure (1) 
we show the structures of three cylindrical single wall 
nanotubes for selected values of m and n. The tt elec- 
tron eigenstates vE'(r) satisfy periodic boundary condi- 
tions on the cylinder, i.e. fy(r + Th) = ^(r). Inter- 
estingly this does not imply periodic boundary condi- 
tions for the eigenfunctions of H c g . Instead, i? c ff is ob- 
tained from a factorization of the single particle state 
"J = ^i(r) • (U a (r), Ub(r)) where (U a , Ub) are the two sub- 
lattice components of the eigenstates of the it electron 
Hamiltonian at the critical K point of the zone, and ip(r) 
is a (slowly varying) eigenstate of H c g . In particular, 
the functions U a and Ub are Bloch functions with crystal 
momentum K = (47r/3a, 0) which do not retain the pe- 
riodicity of the original Bravais lattice, but are invariant 
only under the translations of a y3 x v3 superlattice. 
Since \&(r) is invariant under Bravais lattice vector Th, 
one finds that the function ip is conjugate to U, and ac- 
cumulates a phase e~ lK ' Th under translations Th- 

This phase shifted boundary condition is an awkward 
computational as well as conceptual constraint on the 
low energy theory. However, we observe that this con- 
straint can always be enforced by imposing strictly peri- 
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odic boundary conditions on the wavefunction ip in the 
presence of an effective vector potential a w which satisfies 
2ir J a w ■ dl = —K -Tfi where the line integral is taken on 
a closed orbit around the waist of the cylinder. This vec- 
tor potential can be associated with an elementary flux 
of strength $ = ~K ■ Th/2ir which links the cylinder. 

The azimuthal quantum states which satisfy periodic 
boundary conditions on the surface of the cylinder are 
the cylindrical harmonics e lm ^ with integer m. We adopt 
a coordinate system in which the ^-direction denotes the 
(counterclockwise) tangential direction on the tube's sur- 
face, and the £ direction is aligned along the tube. Then 
for the TO-th channel propagating along the tubule 

2lTV 

H e g = -=-ai{m + $) - iva-id^ (1) 



with spectrum E = ±vyql + (2ir/T h ) 2 (m + $) 2 . Thus, 

for any wrapping where $ is an integer, the accumulated 
phase due to the vector potential can be absorbed into 
the definition of the azimuthal quantum number m. In 
particular, there exists an azimuthal state m = — $ for 
which the "mass" term vanishes, and the electronic spec- 
trum is gapless. This occurs for the wrapped tubules in 
which Th is an element of a regular \/3 x V3 superlattice 
of the original graphite Bravais lattice. For the remain- 
ing two-thirds of the wrapped structures, the minimum 
value of \m + $| = 1/3, so that these structures retain a 
nonvanishing gap AE = Anv /3T^. In a nearest neighbor 
tight binding model, with nearest neighbor hopping am- 
plitude t, one has v = itd/2 so that these primary gaps 
depend inversely on the tube radius R, AE = td/R as 
has already been deduced from numerical work by sev- 
eral groups, ft is noteworthy that the symmetry of the 
Dirac spectrum requires that for the metallic tubules for 
each isospin there is only a single azimuthal branch of 
the electronic spectrum which crosses the Fermi energy, 
independent of the radius of the tubule, and that trans- 
port in a single tube is therefore always governed by a 
single transverse channel. 

In addition to the wrapping constraint described 
above, the local shape of the tubule plays an essential 
role in determining the low energy electronic properties. 
Specifically, curvature and shear in a graphite sheet in- 
troduce variations in the local electronic hopping ampli- 
tudes. Consider fluctuations in the hopping amplitudes 
St a along three nearest neighbor bonds r a at a given site 
on sublattice a. The average value of 5t a simply renor- 
malizes the velocity of the Dirac particle. However, the 
variation from bond to bond introduces a new symmetry 
breaking term into H c s. Inserting St a into the lattice 
Hamiltonian and expanding about K, this term has the 
form SH c g = ( w /2)(a^cr _ + a~a + ), where 

a ± = ij26t a e^ (2) 

a—l 

and — o\ ± i<72- Defining the vector 
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h = J2r a St a /t. (3) 

a=l 

it is straightforward to show af — i(2/3d)(h x ±ih y ). This 
term may be written as an effective curvature derived 
vector potential Bcic — z x h, where z is the local unit 
normal vector. We then have, 

H eS = va ■ (p + a w + a c ) (4) 

A similar expression may be derived for the K' point. 
Equation (4) demonstrates that the effect of fluctuations 
in the bond hopping amplitudes are to displace the sin- 
gular point of the Dirac operator in k space, but not 
to remove it. As an electron propagates on the surface 
of a tube it accumulates a phase from both the winding 
condition (through a w ) and from the local fluctuations of 
the hopping amplitudes (through a c ) which it encounters 
along its path. 

The bond hopping fluctuations can now be deduced 
from the shape of the tubule. Here we consider contri- 
butions from bond length variations, and from the mis- 
orientation of 7r electron orbitals on neighboring sites of 
the tubule. The detailed calculations leading to the re- 
sults given below are straightforward but lengthy j?J , and 
will not be presented here. Instead we focus on the key 
results. One finds that the dependence of the hopping 
amplitudes on bond length for a hop along bond f a can 
be expressed in terms of the metric tensor g^ for the 
curved surface so that 

toalt={P/2#y a rii3ii-tii) (5) 

where (3 = d\nt/dlnd gives the linear dependence of the 
bond hopping operator on bond length. We consider the 
7r orbitals to be oriented along the local normal of the 
tubule surface. On a curved surface, the local normals 
on two neighboring sites are no longer perfectly aligned, 
and this misorientation effect also modulates the hopping 
amplitudes. We find that this effect can be calculated 
using both the metric tensor and the curvature tensor 
Kij [§. The result is 

St a /t = -(l/8y a riK jk K lig kl + ... (6) 

There is an additional contribution which arises from re- 
hybridization of the 7r electron states on the curved man- 
ifold, and it can be derived by studying the effects of 
shape fluctuations on the invariant (n ■ T a )(n' ■ f a ) [Q. 
The essential ingredients for the long wavelength physics 
are already contained in the former two contributions, 
and we now consider them in more detail. 

For a tubule in the form of a right circular cylinder 
we need to specify the metric tensor, curvature tensor, 
and the tipping angle 6 which orients a bond of the hon- 
eycomb network with respect to the % axis" along the 
length of the tube. For the right circular cylinder we have 
gij = Sij and the only nonzero component of the curva- 
ture tensor is = l/R. Thus the pure metric contri- 
bution to a c vanishes, and we only have the orientational 
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contribution. The explicit form may be found using (3) 
and (6) along with the fact that for any three vectors A, 

B, C, EJfa ■ A){f a ■ B)(f a ■ C) = (3/4)Re[A+B+<7+]. 
Expressing a c in the new coordinate system, we find 
a+ = a C £ + ia C Q = —e 3l6 d/16R 2 . Note that this con- 
tribution is proportional to the square of the tubule cur- 
vature, and is unchanged under rotations of the tipping 
angle by 27r/3 as one expects from the symmetry of the 
honeycomb lattice. For the "zig-zag" tubes (as shown 
in figure 1(a)) we have a bond exactly aligned with the 
long axis of the cylinder, so that 9 = and a+ is purely 
real. This means that a c is directed along the circumfer- 
ential direction of the tube. Here the line integral J a c -dl 
around the waist of the tube is nonvanishing (and in gen- 
eral nonintegral) , so that the curvature makes a nonzero 
contribution to the effective mass of the Dirac particle in 
equation (4). For the armchair tubes (as shown in figure 
1(c)) 6 = ir/2 so that the vector potential is purely imag- 
inary, indicating that it is directed exactly along the tube 
direction. For this geometry its only effect to rigidly shift 
the electronic spectrum along the direction in momen- 
tum space. This shift has no physical consequence and 
can be completely eliminated from this geometry by a 
simple gauge transformation. Since <E> is an integer for 
all armchair tubes, and the curvature corrections can be 
removed from the Hamiltonian by a gauge transforma- 
tion, the mass term vanishes and all the straight armchair 
tubules remain metallic. 

For a given cylindrical tubule the total band gap is 
\2v(a w + Re[a+])|, and in figure 2 we display a plot of 
the total gaps predicted for all tubules with radii less 
than 15 A. The plot identifies three distinct families of 
tubules: (a) tubes with primary winding induced gaps 
scaling as 1/R with curvature-derived fluctuations scal- 
ing as 1/R 2 , (b) tubes with vanishing primary gaps, and 
a nonvanishing curvature induced gap (these are shown 
on an expanded scale in the lower panel), and (c) zero gap 
(armchair tubules) for which both the primary gaps and 
curvature induced gaps vanish by symmetry. The data 
predicted within this model provide a strikingly complete 
description of numerical data for these gaps obtained 
from a complete tight binding analysis of these tubes 
employing four basis orbitals per carbon site for each 
of these structures [[To). Thus the Dirac model on the 
curved surface correctly represents the low energy elec- 
tronic physics for this system. We remark that the scales 
of these gaps for tubes of radius sa 10 A are by no means 
negligible and can have an important consequences for 
the low temperature transport properties jll) . 

This model can now be extended to far more complex 
structures which contain fluctuations in the tubule shape. 
Physically, just as the uniform vector potential describes 
a homogeneous mapping of the graphite plane onto the 
tubule surface, a perturbation to the tubule shape pro- 
duces a perturbation in the vector potential which then 
can scatter a quantum particle. Here we will focus only 
the armchair tubes, since these are the only structures 



which are metallic in the absence of shape fluctuations. 
We consider the effects of long wavelength twists and 
bends of the tubule, as shown in figure (3) since these 
are the low energy degrees of freedom for the system. 

We find that even a modest twist can serve as a strong 
scatterer for a propagating 7r electron. The dominant 
effects are introduced through the metric tensor contri- 
butions in equation (5). For a tube subject to a twist 
7 = 8^4) we find a+ = ie 3ie (3(R/2d)j. For the armchair 
tube, a c is directed along the circumferential direction 
and provides a gap in the Dirac spectrum of (3f3'yR/2)t. 
Physically this effect is due to an asymmetric compres- 
sion and dilation of the "axial" bonds on the surface of 
the twisted tubule. Here we find that a twist which ro- 
tates the wrapped graphite structure through an angle 
of 7r over a distance of 1/x introduces a gap of 20 meV at 
the Fermi surface. The contributions from the curvature 
induced misorientation of the n orbitals in equation (10) 
are considerably smaller by a factor d 2 / 12 f3R 2 . 

Interestingly, we find that the coupling to bend is much 
weaker. A deformation with a constant bend but no twist 
does not backscatter a particle. The underlying reason 
for this is that an armchair tube with uniform bend has 
a local mirror plane which preserves the symmetry be- 
tween the "axial" bonds. This ensures that the effective 
vector potential points along the tube, so it is ineffective 
for backscattering a propagating particle. In principle, a 
propagating tt electron can scatter from a spatially vari- 
ation of the bend. However, the first order coupling to 
a spatial derivative of the curvature B^Kqq |l2j vanishes 
due to the azimuthal symmetry of the electronic states at 
the Fermi energy. From the higher order corrections we 
estimate a bend induced gap of order ii? 3 d/A 4 , where A 
is the persistence length of the tube. This is of order 0.01 
meV, which is negligible for the situation of experimental 
interest. 

In a real single wall tubule, one expects that the twist 
can be inhomogeneously distributed along the tubule 
length. (In general, one expects no twist for an isolated 
nonchiral tubule, but the torques on a tubule induced by 
three dimensional packing of these structures undoubt- 
edly induced some degree inhomogeneous twist.) For a 
twisted section of tubule connecting two untwisted arm- 
chair tubes, one can regard the connecting segment as a 
weak link between conducting segments. At sufficiently 
low temperature, backscattering from these defects can 
ultimately lead to localization of a quantum particle. 
However, before this idea can be meaningfully applied 
to these systems, the model will have to be generalized 
to describe the competing effect of interwall quantum 
coherence in three dimensional samples built out of sin- 
gle wall tubes. Nevertheless we remark that recent mea- 
surements on ropes composed of armchair tubes, and on 
mats composed of an ensemble of connected ropes, show 
a resistivity which crosses over from a low temperature 
regime with resistivity decreasing with increasing tem- 
perature to a high temperature regime in which the re- 
sistivity is increasing roughly linearly with temperature 
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p3| . This crossover occurs in the range 10K - 200K (de- 
pending on sample quality and morphology) and since 
these temperatures are in an energy range which can be 
easily accounted for by the shape fluctuations discussed 
above it is tempting to associate this crossover with the 
onset of strong backscattering from quenched disorder in 
the tubule twist. It would be quite interesting to quantify 
this idea by measuring the degree of twist which is actu- 
ally quenched into three dimensional samples composed 
of carbon nanotubes. The transport properties of nan- 
otubes subject to a controlled, but variable, torsion would 
similarly provide important information about these ef- 
fects. Finally, we remark that the twist can be thermally 
excited and provides an important temperature depen- 
dent scattering rate for ir electrons propagating along 
the tubule. 
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FIG. 1. Lattice structures for: (a) the ideal zigzag [12,0] 
tube, (b) the ideal chiral [8,6] tube, and (c) the ideal armchair 
[8,8] tube 
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FIG. 2. Gaps calculated for all right circular tubes with 
radii less than 15 A. The tubes fall into three families: those 
with primary gaps which scale as 1/R (top panel, top curve), 
those with zero primary gap but nonzero curvature induced 
gaps which scale as 1/R 2 (lower curve top panel, and shown 
in the expanded scale in the lower panel) and armchair tubes 
with zero primary gap and zero curvature induced gap. 
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FIG. 3. Deformations of a armchair [5,5] tube. The left 
panel gives the ideal tube structure. The middle and right 
hand panels show the effects of uniform bend and twist on 
the structure. 
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